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ABSTRACT 

A crucial test of the Standard Model is the measurement of electroweak 
gauge-boson scattering. In this paper, we describe a generic parameteriza- 
tion aimed at a realistic simulation of weak-boson scattering at the LHC. The 
parameterization implements resonances of all possible spin and isospin com- 
binations, properly matched to the low-energy effective (chiral) Lagrangian, 
includes leading higher-order effects and contains a minimal unitarization 
scheme. We implement the parameterization in the Monte-Carlo event gen- 
erator WHIZARD and present results for complete partonic cross-section inte- 
gration and event generation. We provide a comparison with the effective 
W approximation that previously has been used for most WW scattering 
studies at hadron colliders. 



1 Introduction 



Exploring the mechanism of electroweak symmetry breaking (EWSB) is the primary focus of the 
upcoming LHC experiments ATLAS and CMS. The simplest explanation, the minimal Standard 
Model (SM), suffers from theoretical deficiencies and does not account for all experimental facts. 
Weakly-coupled extensions of the SM such as its minimal supersymmetric version MSSM are a 
possible solution. All weakly-coupled models contain new particles in the range between about 
100 GeV and 1 TeV that are observable at the LHC. Among them are light scalar states, in 
particular one or more neutral Higgs bosons. 

No Higgs boson has been observed so far, and the LHC will finally decide about its existence. 
If no light Higgs boson exists, we have to consider alternatives to the familiar SM. Models 
without a (light) Higgs boson are strongly coupled, hence much less predictive and more difficult 
to handle theoretically. They need not provide new physics below the TeV region. While 
simple strongly-coupled scenarios such as minimal technicolor tend to be at variance with 
known precision data, more advanced models remain valid, and we are not even close to a 
comprehensive view of the possibilities. 

The theory and phenomenology of strong weak-boson scattering (for reviews, see, e.g., 
Refs. [If2f3P]) has been a subject of active research for more than two decades. Early work on 
a strongly interacting electroweak sector [offi] was motivated by the technicolor paradigm [7] . In 
particular, Bagger et al. [5] considered a collection of benchmark scenarios and their observabil- 
ity at hadron colliders; this study was updated for the LHC parameters in [9|10] . Later work 
focused on the sub- TeV behavior and its extrapolation to higher energies [Illll2|ll3|ll4|ll5|ll6j . 
Studies of WW scattering at lepton colliders are also available [TTlll8f 191120] . More recently, 
interest in WW scattering at the LHC was revived in the context of extra-dimension mod- 
els [2112211231121] . 

Since the LHC will start taking data soon, new and detailed experimental studies are under 
way which prepare for the upcoming analyses at ATLAS and CMS. These have to operate 
on a solid theoretical basis. However, the earlier phenomenological studies mentioned above 
have restricted themselves to particular benchmark models, e.g., the SM, technicolor- inspired 
resonances or specific unitary extrapolations of the low-energy behavior. Non-SM models have 
been treated using simplifying approximations, in particular the effective W approximation 
(EWA) |25j. For the future analysis of real LHC data, it will be crucial to get rid of approx- 
imations and treat the problem with full generality, as far as the physics is accessible to data 
analysis. 

The present paper aims at a practical realization of the strongly-interacting scenario that 
is suited for realistic physics simulation and experimental analysis. To this end, we introduce 
a generic parameterization of weak-boson scattering that includes all resonances allowed by 
spin and isospin with free mass and width parameters. We embed this in the generic effective- 
Lagrangian formalism for electroweak symmetry breaking [26f27] and properly match the high- 
energy region to the low-energy expansion. We include the model-independent part of loop 
corrections to the scattering amplitude. For regulating the high-energy behavior, we adopt a 
straightforward (K-matrix) unitarization scheme. This approach cum grano salis encompasses 
all of the specific models studied earlier. 

The parameterization is extended off-shell in a natural way, and thus can be implemented 
in a parton-level matrix element generator. The SM emerges as a special case. Models can 
thus be studied in the context of cross-section calculation and event generation, and there is no 
need for further approximations. The partonic simulation provides complete six-fermion signals 
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and irreducible background. We have realized this as an extension to the public Monte-Carlo 
simulation package WHIZARD [28.29J, and we present numerical results. 

Beyond partonic cross sections and events, the implementation makes it possible to apply 
parton shower, hadronization, and fast or full detector simulation. This should enable LHC 
analyses of weak-boson scattering to derive solid conclusions from comparing simulation results 
with real data, once the latter are available. 

2 Strong Weak-Boson Scattering 

In this section, we consider a generic no-(light-)Higgs scenario. In the absence of a light scalar 
resonance, weak bosons become strongly interacting in the TeV range (30], and the perturbative 
expansion in the weak couplings g,g' breaks down. To the extent that the corresponding 
scattering processes are observable at the LHC, a measurement of the amplitudes is a probe of 
new physics in electroweak symmetry breaking. 

2.1 The LHC Case 

The LHC can access this kind of physics in processes of the type qq — > jj + 4f. Among the 
Feynman diagrams there are some where the initial quarks radiate approximately on-shell W 
and Z bosons and become hard forward/backward (low-pr) 'spectator' jets, Fig. [TJ The weak 
bosons scatter quasi-elastically and decay into four additional fermions which appear more 
centrally. This is the strong-scattering signal that we are interested in. It depends on detection 
efficiency and background reduction, which W/Z decay modes (four leptons, semileptonic, all 
jets) are useful. 




Figure 1: Topology of vector boson scattering in proton-proton collisions. 

As an alternative to an expansion in the weak couplings g,g', one can expand in powers of 
E/A, where E is the characteristic energy scale of the subprocess, and the cutoff A is loosely 
defined as Attv with the electroweak scale v = {■\/2Gf)~ 1 ^ 2 = 246 GeV. In practice, this 
expansion is valid up to about 1 TeV, where scattering amplitudes approach the saturation 
of unitarity limits. The corresponding effective Lagrangian is known as the electroweak chiral 
Lagrangian [27] . (There is a close analogy with the chiral-Lagrangian approach to low-energy 
QCD [26,31,32].) For each scattering process, the leading order (LO) in this expansion in E/A is 
completely predicted from low-energy data, while the next-to-leading order (NLO) coefficients 
a.i have to be determined by experiment. Some of the parameters have been constrained by 
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Z-pole and W pair production data. LHC data, hopefully, will probe weak-boson scattering 
well into the TeV range and thus provide the information that is still missing. 

A meaningful experimental analysis of a non-perturbative scenario requires a class of models 
to compare with. For each amplitude, the low-energy region which is quantitatively described 
by an effective Lagrangian, has to be matched to the region of unitarity saturation at higher 
energies. In this region, amplitudes may exhibit resonances, or they may approach saturation 
only asymptotically. There is the actual possibility of a rich high-energy structure (like in 
QCD), but we have to keep in mind the limited event rates and energy range of the LHC: while 
the distinction of leading resonances from a structureless amplitude or from each other becomes 
feasible, looking further beyond and determining asymptotic behavior is quite a challenge. 

2.2 Modeling Terra Incognita 

A comprehensive list of phenomenological models for strong EWSB includes all types of res- 
onances that can emerge in quasi-elastic weak boson scattering VV — > VV with V = W, Z. 
The case V = 7 can be ignored: the strong interactions we are interested in are a property of 
the longitudinal degrees of freedom, which are absent for the photon. For similar reasons, we 
do not consider resonance couplings to the other gauge degrees of freedom, i.e. transversally 
polarized W/Z bosons. There is no obvious relation of such effects to electroweak symmetry 
breaking. Similarly, the couplings of a new resonance to SM fermions may be important, but 
with our current knowledge the relation to electroweak symmetry breaking is obscure, so we 
do not take them into account at the present stage. Of course, the model may be extended to 
cover all of these effects as well, if necessary. 

Spin selection rules restrict VV resonances to scalar, vector, and tensor type. In a generic 
approach, resonance masses and widths are arbitrary parameters, with the limiting case M — > 
00 included. For each resonance, the partial width for decay into (longitudinal) vector bosons 
is determined by the couplings to the corresponding scattering channel and sets the lower 
bound for the total width. As stated above, we neglect other couplings, so the VV couplings 
are directly related to the total width. Expanding results for low energies, each resonance 
contributes a calculable shift to the chiral-Lagrangian parameters. 

Low-energy weak interactions approximately respect weak isospin, also known as custodial 
symmetry, SU(2)c [33]. Models with significant violation of weak isospin at high energy tend 
to provide a shift to the low-energy p parameter that is not supported by LEP precision data. 
In this paper, we therefore extend weak isospin to high energies and consider the following 
resonances in VV — > VV processes: 

• scalar singlet a, scalar quintet 0, 

• vector triplet p, 

• tensor singlet /, tensor quintet a, 

with arbitrary masses and widths, including M — > 00. We might also list it (scalar triplet), u 
(vector singlet), etc., but their couplings to weak bosons are isospin- violating and thus either 
small, so we can ignore them, or require unnatural cancellations to preserve the p parameter. 

It is straightforward to classify models of EWSB, also weakly-interacting ones, according 
to their resonance content in VV scattering. For instance, a specific model with a a resonance 
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is the SM. The vector resonance triplet p appears in technicolor models, but also in extra- 
dimension models where it is understood as a W/Z resonance [21J. A tensor / could be a 
graviton resonance [34|. 

We should expect superpositions of resonances. In particular, multiplets with specific 
SU{2) L quantum numbers Ii decompose into superpositions of SU{2) C multiplets: for in- 
stance, the Ii = 1/2 Higgses of the MSSM decompose into a light singlet a = h and a heavy 
triplet 7r = (H + , A, H~). With increasing mass, the latter decouples from VV scattering due 
to isospin. Similarly, the Littlest Higgs model [35] contains a heavy complex 1^ — 1 multiplet 
which decomposes into a scalar 1 = 2 quintet and a singlet. The parameterization that we 
introduce below supports multiple resonances (one per scattering channel). For our numerical 
results, we have switched on only one resonance at a time. 

2.3 Unitarity 

Since we are interested in strongly coupled phenomena in energy ranges where perturbative 
expansions break down, phenomenological models must have unitarity bounds explicitly built 
in. For instance, the LO naive result for the WW — > ZZ on-shell amplitude yields quadratic 
rise with energy, while unitarity at most allows for an asymptotically constant value. In a 
physics simulation, the naive result would produce by far too many events at high energy, while 
in reality there might be no sensitivity to this region at all. 

For quasi-elastic VV —>■ VV scattering, the unitarity requirement is rather simple: the eigen- 
amplitudes, properly normalized, must lie on the Argand circle \a(s) — i/2\ = 1/2. (Strictly 
speaking, this is true in the limit g < £/A where masses are neglected, and photon and 
inelastic channels are considered subleading and are omitted.) For a(s) = 0, this law is trivially 
satisfied. A resonance corresponds to the amplitude crossing the value a(s) = i. 

Conservation of angular momentum implies that the eigenamplitudes have definite angular 
momentum (0, 1,2,.. .), and since the weak bosons have spin 1, at LO there is no unitarity 
problem for angular momentum higher than 2. Furthermore, if we keep weak isospin as a 
symmetry, the eigenamplitudes also have definite isospin quantum numbers. The relevant 
channels coincide with the list of resonances given above. 

Computed at finite order in perturbation theory, a model amplitude that rises from a small 
value of a(s) near s = 0, will eventually depart from the Argand circle. For instance, the 
LO higgsless SM eigenamplitude %q(s) = 2s /v 2 breaks the unitarity limit Re a(s) < 1/2 for 
E > 1.2 TeV, and in a perturbative expansion this is not remedied by loop corrections in finite 
order. Therefore, unitarization models have been invented. They act as an operator that takes 
a scattering amplitude and projects it onto the Argand circle in an ad-hoc way. 

For practical purposes, only gross features of the unitarization scheme are relevant. For 
instance, in ILC physics (y/s < 1 TeV), unitarity saturation is not even reached, so the low- 
energy expansion taken at face value is usually sufficient. The LHC can probe higher energies, 
but both quark and weak-boson effective structure functions fall off rapidly with rising energy 
and strongly suppress the impact of the multi-TeV range. So, the most important property of 
any scheme is that it does ensure unitarity, and thus prohibits any fake s n rise of the amplitude 
that, in a simulation, would produce too many events with large VV invariant masses. 
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3 Basic Theory 



3.1 Effective Lagrangian 

Without a light Higgs boson, the interactions of fermions and vector bosons depend on an 
infinite number of parameters. However, if the S'-matrix is expanded in a series E/A with 
A = 47rv, at any fixed order in the expansion only a finite subset of the parameters is relevant. 
Order by order, the expansion can be generated by a suitable low-energy effective Lagrangian. 

For a useful approximation, the effective Lagrangian has to respect the low-energy sym- 
metries, in particular electromagnetic U(l) and QCD 577(3) gauge invariance, which therefore 
are realized linearly on the fields. The electroweak symmetry SU(2)l X U(1)y is broken by 
fermion and boson masses, but manifest in the low-energy current algebra as well as in the 
massive vector-boson couplings. This can be encoded in a nonlinear realization. Grouping 
quarks and leptons as left-handed and right-handed doublets Ql/r and Ll/ r , one introduces a 
matrix- valued field £(x) which transforms as 

E^U L Xtf R (1) 

under local SU (2) l X U(l)y transformations, where Ul{x) = exp (i Yla=i /? a (x)r a ) and U R {x) = 
exp (if5°(x)T 3 ) with gauge parameters /3 a (x) and Pauli matrices r a . The £ matrix field is also 
a special unitary matrix, i.e., it can be parameterized by 

£(» = exp ^w(x)^ (2) 

with a scalar field triplet w = Yf a =i wa-r< \ cf- App. IA.21 The ground state for the perturbative 
expansion is defined by £ = 1, i.e., w a = 0, and the nonlinearity appears in the w a gauge 
t r ansf or mat ions . 

With these definitions, an effective Lagrangian which generates the lowest order in E/A is 
the chiral Lagrangian [2TB] 

C = ^ tr [(L> M £)t(^£)] - | tr [W^W""] - § tr [B^B M 1 - § tr [G^G""] . 

+ Q L ipQ L + Q R ipQ R + L L ipL L + L R iflL R (3) 



(Q l J:M q Q r + L L EM L L R + h.c.) - L c l E*M Nl ^— £L l - L c r M Nr — - — L 



R 



As the basis for perturbation theory in the gauge couplings g s , g, g', and E/A, this Lagrangian 
accounts for all particle-physics measurements that have been possible so far. 



3.2 Resonances 

To describe resonances in WW scattering, we add new degrees of freedom to the chiral La- 
grangian ([3]): scalar fields o and <fi, a vector field p^, and tensor fields and a^, represented 
by tensor products of Pauli matrices in SU(2) space. In our conventions, they all transform as 
matter fields under SU(2) l according to their isospin representation, 

o^a, p^u LP Ul <t>^(U L ®U L )<t>{U L ®U L )\ (4a) 
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/ and a analogous to a and 0, respectively. 

In terms of physical (charged) fields, the iso-singlets a, f are neutral, 



a 



a 



(5) 



the iso-triplet p decomposes as 



V2 



V2 



P r 



and the iso-quintet fields 0, a contain doubly-charged components, 

(f)=y/2 (0++T++ + 0+T+ + 0°T° + <fTT- + <j>~~T—) 



(6) 



(7) 



where r ++ = r + <g> r+, etc. (App. IA.10 . 

A minimal Lagrangian for these should contain a kinetic term and the lowest order (in a 
derivative expansion) of couplings to W/Z pairs. There are two possibilities: (i) couplings 
to transversal gauge bosons via the field strength W^, B^, and (ii) couplings to longitudi- 
nal gauge bosons via the covariant derivative of the matrix field S. We do not consider the 
first discussed above, such couplings are not directly related to EWSB. Furthermore, 

transversal gauge bosons are associated with a factor g or g' instead of E/A, so the interactions 
of transversal gauge bosons with a resonance are numerically sub dominant. 

Let us look at couplings of a heavy resonance to longitudinal gauge bosons. As shown by 
Appelquist /Longhitano et al. [27], all possible terms can be expressed via the two derived fields 

V M = £( J D / ,£) t and T = Sr 3 S f . (8) 

In the unitarity gauge where S = 1, they reduce to = —igW^ + ig'B^ and T = r 3 . If 
we insist on isospin (custodial symmetry), the isospin-breaking spurion T can be omitted, 
and all couplings to longitudinal gauge bosons proceed via couplings to V M . This vector field 
transforms under SU(2)l x U(l)y as 



Vu^U L V u Ut.. 



(9) 



Each Lagrangian consists of a kinetic term for the resonance and a linear coupling to a bosonic 
current. Explicitly [18] . 



-a {Ml + d 2 )a + aj a 



-tr [cf> {Ml + d 2 ) <t>] +trfofc] 



M 2 i 



M 



^kin ~ fnvf^ ~\~ f^ivjf 



M 2 



(10a) 
(10b) 

(10c) 

(lOd) 
(lOe) 
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(the explicit form of the the M = kinetic term £km of the tensor [App. IA.3j is not needed) 
with the currents 
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. tr [VV] + 



tr [VpV p [ 



(lOf) 

(10g) 
(lOh) 
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tr [V,V] 



The form of the interactions is completely determined by the transformation laws of the fields 
and by the conditions of symmetry and transversality, 



and tracelessness with respect to SU(2) 

tr [pj = tr [0] 



0. 



/Ltl/ 



0. 



<9 M a 



0. 



tr [a, 



0. 



(11a) 



(12) 



Analogous relations hold for the currents and uniquely fix their form, up to terms with higher 
powers of derivatives. 

Higher-derivative terms in the amplitude can be expanded about the resonance location. 
Their on-shell values renormalize the leading interaction terms as given above and can thus be 
dropped. The off-shell corrections are non-resonant and thus renormalize the NLO low-energy 
effective Lagrangian, so they are included there and can also be omitted. In short, our list of 
resonance interactions with longitudinal W/Z bosons is exhaustive (for the vector resonance 
case, see App. [Uj) . 

With the interaction Lagrangian fixed, we can evaluate the partial widths for resonance 
decay into vector bosons. Given the fact the we do not specify couplings to transversal bosons, 
we can only calculate the leading term in the electroweak coupling expansion, which is easily 
computed using the Goldstone-boson equivalence theorem (GBET) The results are 

listed in Table [TJ With increasing number of spin and isospin components, the resonance width 
decreases. Furthermore, with our normalization convention for the dimensionless couplings 
the width of a vector resonance has a scaling behavior different from the others. 
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Table 1: Partial widths for resonance decay into longitudinally polarized vector bosons, computed 
using the GBET. All values have to be multiplied by the factors g 2 /6Air and M 3 /v 2 , where g is 
the coupling in the corresponding resonance Lagrangian. 

In a purely phenomenological approach, the couplings gi in the interaction Lagrangian have 
no meaning on their own, and their normalization is arbitrary. Thus, it is useful to eliminate 
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them in favor of the resonance masses and widths which are observables, using Table [TJ We 
will do this in the following section, so the matching to the low-energy effective theory is made 
free of this ambiguity. 
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= k* 2 </tr [B^[V, 
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3.3 Low-energy effects 

Below the first new resonance, physics is described by the chiral Lagrangian with a double 
perturbative expansion in the electroweak and strong couplings, and in E/A. The LO in E/A 
is generated by the Lagrangian ([3]). The NLO in E/A is generated by one-loop corrections 
and by higher-order operators a,iCi with coefficients a^. The list of NLO terms with isospin 
symmetry SU{2)c consists of [27] 

(13a) 
(13b) 
(13c) 
(13d) 
(13e) 

The first two terms introduce isospin breaking in the same form as the SM, i.e., only via the 
coupling to the hypercharge gauge boson, just as the lowest order Lagrangian does. This 
breaking disappears in the limit g' <C g. C± corresponds to the S parameter, which is well 
constrained by LEP data. £ 2 and £3 affect three-boson couplings and are also constrained by 
LEP; these bounds will be improved by weak-boson pair production at the LHC. The last two 
terms are observable only in weak-boson scattering and are thus unconstrained so far. 

There are several sources that contribute to the a parameters. First of all, they arise as coun- 
terterms for the one-loop correction, and therefore logarithmically depend on a renormalization 
scale. Calculable contributions are generated by integrating out heavy degrees of freedom, in 
particular the resonances introduced above. Ultimately, the values of 014 result from matching 
the underlying theory to the chiral Lagrangian; e.g., in a technicolor model contributions to 
ad can be estimated from technifermion loops. In the analogous case of low-energy QCD, such 
estimates are feasible, while in the electroweak case, the underlying theory is unknown. 

Here, we consider the contributions that result from integrating out resonances at tree level. 
Formally, we can cast the interactions of a resonance $ in the form 



= z 



i$(M 2 + A)$ + $J 
2 



(14) 



with a coefficient z and composite operators A and J. The specific formulae include sums over 
spin and isospin indices. 

Performing the path integral over $, we arrive at the effective Lagrangian which we expand 
in powers of 1/M 2 to obtain 

As far as this Lagrangian contains terms that are already present in the LO chiral Lagrangian, 
they renormalize the LO coefficients, i.e., the couplings g and g' and the electroweak scale 
v. Since the values of these parameters are determined by low-energy data (in the sub-TeV 
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range), those shifts can be ignored. The leading part of the remainder can be expressed as a 
combination of the NLO operators listed above. The resulting contributions to 04 and 0J5 are 
given in Table [2j The values increase with increasing spin and isospin, and expressed in terms 
of the observable parameters v, T and M they all have the same scaling factor v 4 /M 4 . 
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Table 2: Shifts in the NLO chiral Lagrangian coefficients 04 and 05 that result from integrating 
out a heavy resonance at tree level. All values have to be multiplied by the factors I671T jM and 
v 4 /M\ 

If a Lagrangian is used that contains a resonance explicitly, these shifts of the a parameters 
have to be omitted since they are replaced by the low-energy tail of the resonance. Vice 
versa, if the resonance is not explicitly included in the Lagrangian but assumed to be present 
(presumably, because its mass is beyond the reach of the experiment), the oti shifts due to the 
resonance have to be added to the low-energy effective Lagrangian. 




Figure 2: Anomalous couplings 0:4/5 in the low-energy effective theory coming from the different 
resonances under the assumption of identical masses and widths (Tabled). The dashed arrow 
indicates the shift due to renormalization scale variation. (The derivations are given in the 
text.) 

In Fig. [2J we display the directions and relative magnitudes of these shifts in the 0:4-0:5 
plane. We observe that the contributions due to resonances are roughly orthogonal to the shift 
which is attributed to a change of renormalization scale in the one-loop corrections (1301) . which 
makes the two sources distinguishable in principle. Furthermore, arbitrary resonance patterns 
induce a combined shift which lies between the upper and lower-right directions in Fig. [2J This 
coincides with the region favored by causality considerations [To] . 

If there is only one important resonance, a simultaneous fit to both a parameters in the low- 
energy region would thus enable us first to distinguish the isosinglet case (scalar or tensor) on 
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the one hand from the isotrip let /-quintet case (scalar, vector or tensor) on the other hand. If the 
resonance can actually be produced, an angular analysis of its decay products (for instance, in 
the golden channel R — > ZZ — > 4/x) could then distinguish scalar from tensor. The p resonance 
multiplet has the characteristic feature that the ZZ decay channel is absent, a manifestation 
of the Landau- Yang theorem. 

3.4 Reparameterizations 

In this section we discuss alternative parameterizations of the physics we are interested in. Due 
to the equivalence theorem of quantum field theory [57ft. they can lead to different intermediate 
results (such as Feynman rules), but ultimately have to yield the same observables. 

(a) In the previous sections, we have chosen a particular representation of the effective La- 
grangian which manifestly exhibits SU (2) l X U(l)y gauge symmetry and SU (2)c global isospin 
symmetry. While gauging electroweak symmetry is useful for making contact with the SM and 
to low-energy current algebra, and for computing loop corrections, tree-level calculations (at 
least) can be done in unitarity gauge, where weak bosons are merely heavy matter fields. The 
rules for unitarity gauge are 

w^O, (16) 
£ - 1, (17) 
V^-igW^ + ig'B^. (18) 

In this gauge, the Goldstone scalars w a disappear, and only physical degrees of freedom are 
present. 

(b) Alternatively, in the limit that the gauge couplings g,g' can be neglected compared 
with E/A (gaugeless limit), one may omit the gauge fields and study processes with external 
Goldstone scalars w a only. These calculations are particularly simple. Due to the GBET, in 
the gaugeless limit the resulting observables are identical to observables where the Goldstone 
scalars are replaced by physical, longitudinally polarized, vector bosons. 

(c) The UET states that physical observables are invariant with respect to arbitrary non- 
linear field redefinitions. While manifest symmetries should be kept in a linear realization for 
obvious reasons, there is much freedom in the treatment of nonlinear symmetries. A simple 
corollary implies that all parameterizations of the unitary matrix £ in terms of three scalar 
fields are equivalent. For instance, we could alternatively use 

£(*) = ^Jl-^f- x (l - i-wto) (19) 

and get new Feynman rules, but identical results for Goldstone scattering and vector-boson 
scattering observables. 

(d) A straightforward nonlinear reparameterization involves omitting the B field from the 
covariant derivative in (jSJ) and expressing the couplings in terms of 

W£ = tr[V„r a ] and B M , (20) 

1 We may call this the universal equivalence theorem (UET) to distinguish it from the Goldstone-boson 
equivalence theorem (GBET) [36] for electroweak interactions, which is a corollary of the UET and gauge 
invariance. 
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which results in vector fields that are invariant under SU(2)l but transform nontrivially under 
U(1)y instead: W become matter fields while W° behaves like a gauge field. Analogously, by 
multiplying fermion doublets with X factors, fermion fields transforming just under U(l)y can 
be introduced. This approach, which is close to choosing unitarity gauge, has been described 
in Ref .[38]. 

(e) The CCWZ version of the chiral Lagrangian [31] introduces the square root of S, 

£ = ft, (21) 

so in the exponential parameterization 

= exp (^w(z)) . (22) 

The field £ has a mixed transformation law, 

£ - ^4 = UcZUl (23) 

which defines an 577(2) matrix function of the transformations Ul(x) and Ur(x) 

and of the field £(x). The matrix Uq{x) can be interpreted as a local isospin transformation, 
C/ c e 5[/(2) c . 

Using £, the chiral fermion multiplets Ql/_r and Ll/_r can be promoted to Dirac spinor 
multiplets, 

which no longer transform under SU(2)l or f/(l)y, but have a common transformation law as 
isospin doublets: Q — > £/cQ, £ - ► t^c-^- Similarly, £ factors make the resonance multiplets 
invariant under SU(2)l x C/(l)y, but transforming under SU{2)c- 

For a vector resonance p, the CCWZ formulation allows to introduce it either as a matter 
field, or as the gauge field of local SU(2) C , with gauge couplings only. In the Lagrangian above, 
we have introduced the p resonance as a matter field. In App. (UJ we describe the alternative 
formulation with p as a gauge field and verify the equivalence of the two approaches. 

To summarize, while our formulation of the chiral Lagrangian coupled to resonances is by 
no means unique, it is nevertheless equivalent to any other formulation that correctly describes 
low-energy physics. As such, the chiral Lagrangian approach is model-independent. We do 
use model assumptions and truncations, however: no isospin violation beyond hypercharge and 
fermion couplings, minimality in the number of degrees of freedom (at most one resonance per 
channel), a minimal set of couplings (no independent couplings to transversal gauge bosons, 
no self-couplings of resonances), truncation of the low-energy expansion (LO and NLO only), 
and minimality in the unitarization scheme (no extra parameters) . As long as the new degrees 
of freedom are heavy, these model assumptions are likely irrelevant for the experimental pre- 
cision that can be achieved at the LHC. Extensions of our approach, e.g., including secondary 
resonances, are easily possible, but not worked out here to keep this paper compact. 

In App. |D] we relate various specific models that are frequently used in the analysis of 
weak-boson scattering to our generic parameterization. 



11 



4 On-Shell Scattering Amplitudes 



4.1 Low- energy effective theory 

Let us look first at the W + W~ — > ZZ weak-boson scattering amplitude. In the electroweak 
coupling expansion, the leading term is of order g° and corresponds, at high energy, to the 
scattering of longitudinally polarized particles. This term rises with s, while the scattering 
amplitudes of transversally polarized vector bosons come with factors of g and asymptotically 
do not rise with energy. By the GBET, the leading term is equal to the amplitude A(s, t, u) for 
w + w~ — > zz Goldstone scattering. This amplitude is easily computed using the Lagrangian (j3J). 
At tree-level, but to NLO in the E/A expansion, it is 



.4 



tree ( 



S, t, u) 



+ 4<x 



t 2 + u' 



+ 8a 5 - 



V* V* V* 

The leading real part (order g°) of the one-loop correction is given by [39] 



(25) 



.4 



l-loop 
C 



(s,t,u) 



1 



16tt j 



1 



- In fr + 8C 5 — + 



t(s + 2t) /! 



6v 4 



ln^ + 4C 4 - )+(t 



t 



u 



(26) 



where fi is the renormalization scale, and C 4 and C5 are finite scheme-dependent matching 
coefficients. For instance, in the MS scheme, /1 is identified with the MS scale, and C4 = C5 = 0. 
By contrast, in the scheme where a fictitious (heavy) Higgs boson is used as a regulator [40J, 
we have 



M 



11 



and Ca 



1 

18 



-0.056, C 5 



9tt 
16^ 



37 
36 



-0.0075. 



(27) 



Note that these matching coefficients are numerically small, so the difference between the two 
schemes may be neglected. Other schemes are possible, e.g., the QCD-inspired scheme used in 
Ref. PS] is reproduced by C 4 = -13/72, C 5 = -5/72. 

In Fig. [3J we plot the angular dependence of the one-loop correction. If the renormalization 
scale fi is chosen equal to the energy y/s, the loop correction, and thus the angular dependence, 
is less than 2.5 %. Since the NLO correction is proportional to s 2 (compared with the LO 
amplitude proportional to s), it rapidly becomes important for s > /1 2 . However, this mainly 
indicates the breakdown of the low-energy expansion at high energies. 

We can transfer the scheme-dependent matching coefficients to the NLO counterterms, so 
the above result is reproduced by maintaining only the logarithmic terms in the amplitude, 



16tt 5 



2v 4 



In 



M 2 t{s + 2t) . M 



+ 



Qv 4 



In 



u 



and adding one- loop matching contributions to a 4 and a 5 , 



a 



(i) 



1 



2 C 4 , 



a: 



(i) 



16tt 2 ~" "° 167T 2 

The renormalization scale dependence of these coefficients is given by 



a. 



a 4 (/,) = a 4 M-^ I ^ln^, 



asO) = as(^o) - ^7 In ^ • 
24 167r^ /1q 



(28) 



(29) 



(30) 
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Figure 3: Dependence of the loop correction on the scattering angle 9. The parameters are 
^ = fj, = 1 TeV. 



with some reference scale fiQ. 

Isospin symmetry determines all individual scattering amplitudes in terms of the master 
amplitude A(s, t,u): 



A{w + w~ —* zz) = A(s,t,u) 

A{w + z — > w + z) = A(t, s, u) 
A{w + w~ — > w + w~) = A(s, t, u) + A(t, s, u) 
A(w + w + w + w + ) = A(t, s, u) + A(u, s, t) 

A(zz — > zz) = A(s, t, u) + A(t, s, u) + A(u, s, t) 



(31a) 
(31b) 
(31c) 
(31d) 
(31e) 



Expanding the amplitudes in powers of the energy, the order- E 2 term is known as the low-energy 
theorem (LET) gT]: 



A (0) (w + «r -> zz 
A {0 \w + z -> w + z 
A^°\w + w^ — > w + w~ 
A {0 \w + w + w + w + 
A {0 \zz -> zz 



s/v 2 

t/v 2 

—u/v 2 

-s/v 2 





(32a) 
(32b) 
(32c) 
(32d) 
(32e) 



These expressions are mo del- independent and depend just on the electroweak scale v. 



4.2 Resonances 



In Sec. 13.21 we have introduced heavy resonances in weak-boson scattering. The interaction 
Lagrangians (jlOaHlOel) induce couplings to vector bosons and to Goldstone bosons, which are 
related by electroweak gauge invariance, maintaining the GBET. Each resonance multiplet 
therefore contributes additional terms to the Goldstone scattering amplitude A(s,t,u), which 
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have poles at the appropriate locations. We do not yet include the resonance widths. The new 
contributions are 

^( S>i ,„) = _|_^_ (33a) 

+ ^ + (330 

A '(^ = -i^ P ^ + W (33d) 
a 2 f t 2 u 2 /2 s 2 s 2 \ 1 

"•(«. ») = { tnss^ s - "> + ^ P2(u - *■ () - (st^f - p ^ «>} 

(33e) 

where P 2 (s,t,u) = [3(t 2 + u 2 ) - 2s 2 }/ s 2 . 

Beyond the resonance location, for g a — 1 the a exchange amplitude cancels the rise of the 
LET amplitude. This is the SM case. Otherwise, beyond the resonance all amplitudes rise with 
a power of s/M 2 . This implies again unitarity violation, which has to be cured by the unknown 
UV completion of the theory. 

4.3 Eigenamplitudes 

For the analysis of unitarity, we need the spin-isospin eigenamplitudes, i.e., scattering ampli- 
tudes for superpositions of states which scatter only into themselves. We first list the isospin 
eigenamplitudes 

Aq(s, t, u) = 3A(s, t, u) + A(t, s, u) + A(u, s, t) (34a) 

A ± (s, t, u) = A(t, s, u) - A(u, s, t) (34b) 

A 2 (s, t, u) = A(t, s, u) + A(u, s, t) (34c) 

which can be decomposed into partial waves using Legendre polynomials, 

oo 

Ai(s, t,u) = J2 A u( s ) (2 ^ + 1) Pj(s, t, u), (35) 
j=o 

where Ajj ^ only for I — J even. The coefficient functions Ajj(s) are the spin-isospin 
eigenamplitudes. They are obtained by angular integration, 

r° dt 

Au(s)= / -A I (s,t,u)P J (s,t,u). (36) 

J -8 S 

Below, we explicitly list the spin-isospin eigenamplitudes, treating LO, NLO, and resonances 
separately: 

(a) The eigenamplitudes for the LO Lagrangian: 

a(°)-o s a {0) - S 4 (0) - S (17) 

A 00 -*~2 A ll -3^2 A 20 - ~2 y 6 ') 
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All other terms vanish at this order. 

(b) The one-loop correction with its logarithmic angular dependence contains partial waves 
of arbitrary spin. We extract the leading logarithms ln(/z 2 /s), project out the partial waves 
and truncate the series at spin 3, which numerically is an excellent approximation. Adding the 
tree-level NLO coefficients, which should include their scheme-dependent and scale-dependent 
parts ([251 S5D , the real part of the result is 



-4, 



(i) 
oo 



^02 



A (l) 



g 

- (7a 4 (» + lla 5 (/i)) + 

8 1 

— (2a 4 (fi) + aM) + — 



'4 1 

- (a 4 (n) - 2a s (fj,)) + 



+ 



7 



16tt 2 V 1080 



.4 
.4. 



(i) 

20 

(1) 
22 



16 1 

— 2a 4 /i + a 5 /i + -— 

' 4 1 

— (a 4 (^) + 2a 5 (^)) + Y^2 



25, u 2 11 
— In — H 

9 s 54 

9 s 135 



1 

54 



10, u 2 25 
— In — H 

9 s 108 



45 s 



247 \ 



5400 



(38a) 
(38b) 
(38c) 
(38d) 
(38e) 
(38f) 



We note that the scale dependence of the a parameters cancels the scale-dependence of the one- 
loop terms, as it should be the case. The results are shown in Fig.Hl While the loop correction is 
small below about 1 TeV, for higher energies it becomes important and, eventually, drastically 
changes the behavior. For instance, in A 0Q there is a cancellation between the LO and NLO 
terms at 2 TeV. This clearly indicates the breakdown of the low-energy expansion. 




[GeV] 



Figure 4: Spin-isospin eigenamplitudes as functions of the energy yfs, unitarized according to 
the scheme described in Sec. \4-5\ Solid curves: LO; dashed curves: NLO. Dotted curve: A 00 (s) 
without unitarization. 
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(c) For the decomposition of resonance corrections, we define the following functions: 

r° dt t 2 

Sj(s) = j -—^p (jt, s ,u)Pj(s,t,u) (39a) 
f° dt t 

Vj(s) = J -j-j^P 1 {t,s,u)P J {s,t,u) (39b) 

r° dt t 2 

Vj{s)= j—^P 2 {t,s,u)Pj{s,t,u) (39c) 

*/ — S 

which we give explicitly in Appendix IA.4[ We obtain for the isosinglet scalar, 



^00 1 


» = 


9l s 2 
v 2 s-M 2 




» = 


Q 2 

-2%S 2 {s) 


A° n \ 


» = 


V 



- 2-M s ) 

v 



AUs) = 


-2%«S 3 ( S ) 


(40a) 


A° w {s) = 


q 2 


(40b) 


A 22 (s) = 


-2^S 2 (s) 


(40c) 



the isoquintet scalar, 



59i 



<o(s) = -^So(s) Af 3 (s) = -^S 3 (s) (41a) 

Sv z or 

2 2 2 2 

<M = - 5 -%S 2 (s) AUs) = Ah-^TH ~ \~Ms) (41b) 

2 v 1 s — M l 6 ir 

<(s) = |§ AUs) = -\^S 2 {s) (41c) 
the isotriplet vector, 



A p 00 (s) = -4 9 2 p V (s) - 3g 2 p ^- 2 A^s) = -2gf-^^S,{s) (42a) 

9 c 4- M 2 s 

AM = -H —^r- «.(») ^SoM = 2 9 2 P„( S ) + (42b) 

« = "f - <£jj}J - 2 ^.M = tfF^SA') (42c) 



the isosinglet tensor, 



2 2 2 

= -^(.) - ^ (43a) 

^w = -^w+sp^ (43c) 

4>(») = "4 fl + 6^ + S 3 (s) (43d) 



3t> 2 V M 2 M\ 
- ^ f 1 + + «4) &(.) - (43f) 
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and the isoquintet tensor, 



AO. 



Aa 
^02 



^20 ( S 



^22 ( S 



5 gl - 
"a ^n, 

OZT 

,2 



108 v 2 M 2 

-.2 \ 



3a S 2 



216 M 2 



9l 



12 3v 2 
5 



^i(s)- 



12 3v 2 
1 



432 w 2 M 2 

2 



St 



12 3u ; 



X>o(s) - 



gl 



10f 

i 



w 2 M 2 



<? 2 



60w 2 s — M 2 12 3v 2 



1 + 6 ^ + 6 ^)*<*> 



2160 v 2 M 2 



(44a) 
(44b) 

(44c) 
(44d) 

(44e) 
(44f) 



The coefficient functions Ajj contain poles in s — M 2 as well as finite parts. The poles are con- 
fined to those (/, J) combinations which correspond to the (I, J) assignments of the resonances. 
Again, we truncate the partial- wave expansion at J = 3, so for each spin-isospin combination 
we only keep the leading and one subleading term. 

4.4 Unitarization scheme 

Elastic unitarity requires that the normalized eigenamplitudes 



respect the Argand-circle condition 



which can also be stated as 



\a IJ {s)-i/2\ = 1/2, 
1 



Im 



au(s) 



-1. 



(45) 



(46) 



(47) 



Computed in finite-order perturbation theory, or deduced from some model, the amplitude a(s) 
will usually fail this requirement. However, an arbitrary amplitude a(s) can be transformed 
into a unitary amplitude if we take the real part of l/a(s) and add —i as the imaginary part, 
i.e., 



a{s) = 



1 



Re(l/a(s)) - % 
a(s) 



if a(s) is real. 



1 — ia(s) 

For the unnormalized eigenamplitudes A IJ (s), this can be rephrased as 



(48) 
(49) 



A u (s) = Au(s) + AA u (s), where AA u (s) = 



AiAs) 2 



327r 1 " 3k A u( s ) 



(50) 
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a(s) 



Figure 5: K matrix construction for projecting a real scattering amplitude onto the Argand 
circle 

This is the K-matrix unitarization scheme [12] , cf . Fig. 

With this prescription, a LET amplitude A(s) = s/v 2 becomes 

A(s) = ™ 32™, (51) 

so instead of rising quadratically with energy, the absolute value of A(s) asymptotically ap- 
proaches saturation, formally a resonance at infinity. 

The K-matrix scheme transforms a simple-pole amplitude, A(s) = —c/(s — M 2 ), into Breit- 
Wigner form, 

A(s) = ^ with T = — — , (52) 

K J s- M- ■ i.\IV 32ttM' v ; 

so it is an alternate implementation of Dyson resummation for s-channel particle exchange. If 
c is not a constant but depends on energy, we get a Breit-Wigner resonance with s-dependent 
width. In particular, the amplitude 

A ( s ) = o ^7 53 

y 1 v 2 s-M 2 y ' 



is transformed into 



M 2 s M 2 

Ms) = --T Z I 777^ s with T = -—M. (54) 



v 2 s — M 2 + jMr^p 32tto 2 
Eq. ([331 has the low-energy expansion 

An expansion of this form can also be treated by the inverse-amplitude method (IAM) for 
unitarization [13]. The result is 

i(s) = AW(s)-AW(s)-^AW(s) 2, (56) 

which equals the (1,1) Pade approximant, and precisely coincides with (1541) . We observe that, 
in the present context, the IAM or Pade unitarization scheme is a special case of the K-matrix 
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scheme, where the low-energy expansion of the amplitude is identified with the low-energy tail 
of a single resonance. In QCD, where the p meson dominates form factors at low energy, this 
turns out to be a valid assumption which leads to accurate high-energy extrapolations. In 
the electroweak case, physics may be different, and the actual (unitary) weak-boson scattering 
amplitudes need not follow the extrapolation of the K-matrix/IAM/Pade or any other given 
unitarization scheme. 

In QCD, low-energy parameters can be computed, to good accuracy, by integrating out the 
p resonance. This may also be the case for the leading resonances in electroweak interactions 
(we list the necessary formulas in Sec. 13. 31) . but there may well be extra contributions that can 
be assigned to further resonances, or to other physical effects. For this reason, we keep «4 and 
«5 as independent parameters in our implementation. 

The detailed shape of resonances in weak-boson scattering may also differ from the (running- 
width) Breit-Wigner that our parameterization provides. However, the experimental resolution 
of weak-boson pair invariant masses at the LHC will be limited, so there is little hope for 
precise resonance scans. A parameterization in terms of mass and width, augmented by extra 
«4 5 5 parameters which describe deviations in the low-energy tail, is sufficient. 

Beyond a resonance peak, our expressions suggest a definite prediction, such as a new rise 
of the amplitude with a definite power of s. We should emphasize that this is misleading: the 
behavior in this region is arbitrary and can only be modeled, introducing further parameters. 
However, any precise measurements of the high-energy tail of a heavy resonance will be chal- 
lenging, if not impossible at the LHC. The only property of unitarized amplitudes that we 
really make use of is: that they do not exceed the unitarity limits. 



4.5 Unitarized Amplitudes 

In this section, we apply the unitarization scheme defined above to the generic parameterization 
of scattering amplitudes. Collecting everything, each eigenamplitude consists of a LO (LET) 
part, a NLO correction which includes the one- loop part and finite extra contributions to the 
a parameters, and resonance terms: 

Aij(s)=Af](s)+4)(s)+ Yl A ?J& ( 57 ) 

R=(7,cj>,pJ,a 

which we write in the form 

Au(s) = Afj(s) + Fjj(s) + (58) 

where Fu(s) is finite, and Gu(s) is proportional to s (vector), or s 2 (scalar, tensor). According 
to the prescription in the previous section, the unitarized amplitude becomes 

Ajj(s) = A \ J ^ = Afj{s) + AAjj(s), (59) 

where the correction to the LET amplitude is given by 

(i s — A/f 2 \ 

1 + ^ m ''> + AG „(,)-(,-^)[l- A (A W (.) + F„(,))] j (60) 
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In Fig. [6] we draw the absolute values of the resulting unitarized eigenamplitudes, including 
the LET part Afj ([37]). Since the resonances have definite spin and isospin quantum number 
assignments, each plot contains exactly one curve with a resonance, while the other curves are 
non-resonant. The resonance masses Mr (R = a,(f), p, f,a) have been set to 1 TeV, and the 
couplings gpi to unity. The unitarization prescription smoothly cuts off the amplitudes, so their 
absolute values do not exceed the limit 32n « 100. Some of the amplitudes (e.g., Aq ) contain 
terms rising like a power with the energy and eventually saturate this bound, while others (e.g., 
Al 3 ) rise logarithmically at most, so at accessible energies they stay much below this limit. 

For the scalar isosinglet a, the choice g a = 1 corresponds to the SM with a heavy Higgs. 
In this case, unitarity is restored already by the scalar resonance exchange. Hence, as long as 
M a is below about 1.2 TeV, the asymptotic values of all A u u stay below the limit of 32n. At 
tree level, they are constants that depend on the ratio M 2 /v 2 . This is slightly modified by loop 
corrections and by the unitarization prescription. For g a ^ 1, the cancellations are incomplete, 
and the amplitudes A a u behave in the same way as the other amplitudes. 

Several of the curves exhibit a zero, which in the logarithmic plots manifests itself as a sharp 
down-pointing spike. In fact, in our parameterization this happens for all resonant amplitudes, 
with the exception of the SM Higgs case. The reason is negative interference between the 
resonant propagator and the contact term; the latter is necessary for satisfying the LET and 
rises with a higher power of the energy. For vector resonances, cancellation typically occurs 
at very high energies (above 10 TeV), while for tensor resonances the effect is visible in the 
energy range that we have chosen for our plots. However, if such a zero occurs beyond the 
resonance mass, it should not be taken seriously, because in this range the amplitude contains 
further, undetermined contributions, and the energy behavior of the contact term as given by 
our formulae is not a prediction. Only if this zero appears below the resonance a dip should 
actually be expected. This is the case for A 2Q in the presence of a scalar isoquintet. 

The analytic behavior of the amplitudes is transparent if we plot the real part, which vanishes 
on a resonance. This is illustrated in Fig. [7J All curves cross zero at 1 TeV, the resonance 
mass. Beyond this, they rise and asymptotically approach zero again. This is the resonance 
at infinity generated by the unitarization procedure. The exception is the a resonance which 
approaches a constant, since in this model (the SM) there is no unitarity problem. 

For a concrete Monte-Carlo implementation, we need the unitarized amplitudes for physical 
states, e.g., w + w~ , zz, etc. Therefore, we first translate the spin-isospin eigenamplitudes back 
into corrections to the isospin eigenamplitudes as functions of s, t, u, 

AA (s, t, u) = AA 00 (s) P (s, t, u) + AA 02 {s) 5P 2 (s, t, u), (61a) 
AAi(s,*,it) = AA 11 (s)3P 1 (s,t,u) + AA 13 (s)7P 3 (s,t,u), (61b) 
AA 2 (s, t, u) = AA 20 (s) P (s, t, u) + AA 22 (s) 5P 2 (s, t, u), (61c) 

The result is shown in Fig. [81 The plot clearly exhibits the characteristic angular dependence 
of the resonances with J = 0, 1, 2, respectively, while the continuum background that we have 
included is negligible for s = M 2 . The nonresonant part is important, however, to describe the 
off-peak amplitude behavior. 
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Figure 6: Unitarized spin-isospin-eigenamplitudes for Goldstone-boson scattering. In each plot, 
we display the eigenamplitudes for a definite spin and isospin value, one curve for each of 
the five possible resonances a,4>, p, f,a. The resonance masses are fixed at 1 TeV ; and their 
couplings to Goldstone bosons have been set to unity. 
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Figure 7: Real part of the eigenamplitudes \Au(s)\, each with the corresponding resonance(s) 
switched on; Mr = 1 TeV. 




12 3 



9 [rad] 



Figure 8: Angular dependence of the amplitudes \Aj(s,t,u)\ for I = 0,1,2, each with the 
corresponding resonance(s) switched on and evaluated at y/s equal to the resonance mass. 
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This, in turn, is translated into corrections to the individual scattering amplitudes, 

AA(w + w~ — > zz 

AA(w + z — > w + z) 

AA(w + w~ — > w + w~ 
AA(w + w + — > w + w~ 
AA(zz — > zz) 

Unit arizat ion breaks crossing symmetry, since it is applied only in the s-channel. Explicitly, we 
obtain 

r «A 

AA(w + w~ — > zz) = 





[s,t,u) 


— 


\aa 2 


's,t,u) 


(62a) 




[s,t,u) 


+ 


\aa 2 


's,t,u) 


(62b) 


\ AA ° 


[s,t,u) 


+ 




s,t,u) + ^AA 2 (s,t,u) 




(62c) 


AA 3 (s 


,t,u) 






(62d) 


l AA ° 


[s,t,u) 


+ 


2 A A 

3 


[s,t,u) 


(62e) 



a 5 + — (AA)o(s) - AA 20 (s)) - — - (AA 02 (s) - AA 22 (s)) 



+ 4 

AA(w + z -> w + z) = 4 



24s 2 
5f 4 

a 4 + ^ (AA 02 - AA 22 ) 



12s 2 



t 2 + u 2 



v 1 
(63a) 



«4 + ^AA 20 (s)-— AA 22 (s) 



AA(w + w~ 



+ 8 
+ 4 
w + w~) = 4 



3f 4 • i , . 15f 4 . . . . 

3t> 4 15f 4 
«4 + ^AAn(s) + — A4 2 (s) 
8s 2 8s 2 



(63b) 



« 4 + 2a 5 + ^ (2AA 00 (s) + AA 20 (s)) 



5v 4 



12s 2 



(2AA 02 (s) + AA 22 (s)) 



+ 4 
+ 8 



a A + 2a 5 + —, (10AA 02 (s) - 3AA n (s) + 5AA 22 (s)) 



8s 2 



aA+ 16^ ( 10AA ° 2 ^ + 3AA n( s ) + 5A ^2 2 (s)) 



u 



(63c) 



W W 



a A + —; {AA 20 {s) - 10AA 22 (s)) 



+ 4 



8s 2 

15f 4 

a 4 + 2a 5 + ^2-AA 22 (s) 



t 2 + m 2 



(63d) 



AA(zz -> zz) 



« 4 + a 5 + 24^2 ( AA oo(s) + 2AA 20 (s)) 



(AA 02 (s) + 2AA 22 (s)) 



12s 2 

5f 4 

«4 + « 5 + g^J (^^02(S) + 2AA 22 (S)) 



t 2 + M 2 



(63e) 



Here, the coefficients functions AAjj(s) are determined by decomposing the results from Sec. 14.31 
according to (1591) and inserting this into the unitarization formula (1601) . 
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Adding the above correction terms to the LET scattering amplitudes f)31aH31e| . we have a 
complete and unitary parameterization of on-shell Goldstone scattering. The parameterization 
depends on and a^, on a renormalization scale fi, and on the mass and width parameters of 
the five possible resonances. 



4.6 Off-shell Implementation 

For realistic calculations, we want to transform the unitarized Goldstone scattering amplitudes 
into matrix elements for off-shell weak-boson scattering. We first note that the complete SM 
without a Higgs and without anomalous couplings, already yields the LET result for weak- 
boson scattering. In order to avoid double-counting, we therefore have to remove the LET part 
from any extra contributions that we add to the theory. This is achieved by considering only 
the correction terms f)62aH62ej) instead of the complete unitarized amplitudes. 

The chiral Lagrangian with NLO parameters (i.e., and 0:5) provides an off-shell for- 
mulation for the low-energy effective theory. We can determine Feynman rules and compute 
complete matrix elements of 2 — > 6 fermion processes which include weak-boson interactions 
with anomalous couplings. The Feynman rules of four-boson couplings depend on a± and a$. 
In unitarity gauge, they are derived from the quartic gauge interactions 



£qgc = e 2 [gTA»A»W-W+ - g^A^W^W+] 

gfA^Z" {W;W+ + W+W-) - 2g?A»Z lt W-"W; 



2 *-"w 



e 2 ° y 
r 2 



e 

+ 



e 2 



+ e 2 ^ [gf z Z"Z v W-W+ - g zz Z^Z ll W~ u W+] 

, gV W W^W^W-W+ - g™ w (W^W;) 2 ] + -^h zz {Z»Z,)\ (64) 

where the SM values of the couplings^) are given by 

g vv' = gVV' = l (yV = 77 , 7 Z , Z Z \WW) , h zz = 0. (65) 

If we include the dependence on all five isospin-symmetric NLO chiral parameters ctj fll3aH13e| . 
the deviations from the SM values are 

Ag{< = Ag? = Agf = Agf = -/^a, + ^a 3 (66a) 

Ag zz = 2Agf + 4«4 Ag zz = 2Agf - ^a 5 (66b) 

AgY W = 2clAgf + g 2 a 4 Agf w = 2^Agf - g 2 (a 4 + 2« 5 ) (66c) 

h zz = g 2 (a 4 + a 5 ) . (66d) 

We can now construct a generic off-shell parameterization of weak-boson scattering that 
corresponds to the unitary on-shell Goldstone scattering amplitudes f l62aH62e1) . When the 
Feynman rule for a given quartic gauge vertex is inserted in a physical process, we replace the 



2 In these expressions, the numerical values of the electroweak gauge couplings and the weak mixing angle 
depend on the precise definition of the electroweak renormalization scheme. 
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dependence on the constant parameters «4 and by form factors which depend on s. For 
W + W~ — > and W + W + — > W + W + , this is straightforward: the two expressions s 2 /t> 4 and 
(t 2 + u 2 )/v A correspond to the two different interaction terms for WWZZ (WWWW) in the 
Lagrangian (I64p . so we simply can replace the constants «4 and (0:4 and 04 + 205) by the full 
expressions ( I62a[ I62d[) . respectively. For the other processes, this assignment cannot be done 
in the interaction Lagrangian, but it is obvious in the Feynman rule, where each term s 2 /v A , 
t 2 /f 4 , and u 2 /v A corresponds to a definite combination of g^ u Lorentz factors. 

As a first result, we can compute on-shell scattering amplitudes for physical W and Z 
bosons. These combine the features of the chosen resonance model with SM effects such as 
photon and W/Z exchange. Since on-shell initial vector bosons cannot be prepared in practice, 
we defer this discussion to Appendix [El 

Such an algorithm breaks crossing symmetry, but this is natural since the unitarization 
scheme already breaks crossing symmetry. In a practical implementation, for a given vertex 
we implement all possible orientations of the time arrow as alternatives, and determine the 
orientation that is actually realized when we insert the vertex into a physical process. This is 
straightforward to do for an automatic matrix-element generator. 

Two sources for ambiguities appear in this construction, (i) The GBET relates Goldstone 
scattering amplitudes to weak-boson scattering amplitudes only in the high-energy limit, and 
only for longitudinal polarization. We do not specify couplings to transversal gauge bosons, 
which are not directly related to EWSB and formally subleading in the physics of strongly 
interacting weak bosons. Corrections to the GBET therefore can be computed only up to 
further free parameters. Keeping this in mind, we translate the Goldstone amplitudes to 
weak-boson amplitudes using the leading-order GBETjf] (ii) Strictly speaking, the Mandelstam 
variables s, t, u in form factors are defined for on-shell scattering of massless particles, t, u can 
be replaced by Lorentz factors which are unambiguous, but in the off-shell continuation, the 
subenergy squared s is evaluated for massive off-shell W j Z bosons. This affects the unitarization 
corrections, but these are scheme- dependent anyway. Their main property - to cancel any 
unphysical rise of subamplitudes - is preserved off-shell. It also affects the location of resonance 
poles. However, as discussed in Sec. 13.21 off-shell effects in the latter are accounted for by 
higher-dimensional operators and translate into corrections to 0^5. Finally, we recall that the 
off-shell continuation of W/Z propagators is controlled by electroweak gauge invariance. We 
keep SU{2) L x U(l)y symmetry manifest in the gauge and fermion sectors by using covariant 
derivatives, so this is not an issue. 

As a cross-check, we can compute 2 — ^ 6 fermion processes for the ordinary SM with a 
Higgs boson. In our parameterization, this is the chiral Lagrangian with 04 = 05 = and a 
a resonance with g a — 1. The form factors for the WWZZ, WWWW , and ZZZZ vertices 
contain exactly the Higgs propagator factors that we would have obtained with the Higgs boson 
as an ordinary particle. In the s-channel, the propagator pole turns out to be regularized by a 
running width F6(s) x s/M 2 , which is a sensible treatment of the width of a heavy Higgs boson 
in SM scattering amplitudes [44J. So, despite the fact that we have used the leading-order 
GBET, our off-shell formulation exactly reproduces the tree-level SM result, both on-shell and 
off-shell. The only missing parts are double-Higgs and Higgs-fermion couplings (see e.g. |45j). 
but those couplings do not contribute to the processes we are interested in. 

3 If we adopt the EWA, part of the difference w.r.t. complete amplitudes is formally of the same order as 
the GBET corrections. However, the EWA strongly affects kinematics and ignores a large set of irreducible 
background diagrams, so the numerical impact of this approximation for LHC analyses is much more important 
than model-dependent ambiguities in corrections to the GBET. 
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5 LHC Processes 



5.1 Monte-Carlo simulation 

We have implemented our parameterization of vector-boson scattering in the multi-particle 
event generator WHIZARD [28,29]. The program generates matrix elements for partonic pro- 
cesses via optimized helicity amplitudes while avoiding the redundancies inherent in a Feynman 
diagram expansion. These optimized matrix elements together with a highly efficient phase- 
space setup enable the simulation of six and eight-particle final states. WHIZARD contains the 
infrastructure for simulations of complex collider environments like structured beams, parton 
shower, and interfaces to fragmentation and hadronization. 

As the starting point for the implementation in WHIZARD, we have chose the SM extension 
with anomalous three-boson and four-boson couplings which has been used for the simulation 
of anomalous triple and quartic gauge operators [20|46II47| . The algorithm for the symbolic gen- 
eration of the matrix elements in WHIZARD, which is especially suited for the inclusion of beyond 
the Standard Model (BSM) physics |48| . allows for the insertion of operators in specific time 
directions necessary by the crossing-symmetry breaking effects of the K-matrix unitarization 
prescription. 

5.2 Comparison with the Effective W approximation (EWA) 

In 2 — > 6 fermion processes that contain weak-boson scattering (Fig. [T]) the W/Z bosons that 
initiate the interaction are represented by their propagators with a spacelike momentum. The 
main contribution comes from the region with small virtuality, and we are interested in the 
region of large cm. energy of the vector boson pair. In this region, the virtualities and the 
masses of the vector bosons induce only small corrections to the amplitude, so the initial vector 
bosons can be treated as approximately on-shell. 

We can thus approximate the dominant Feynman graphs by a convolution of massless split- 
ting (of the initial quark into a quark and a vector boson) with the vector-boson interaction, 
which is called effective W approximation (EWA) |25j : 

o-(qiq2 -> q'tq'XVl) « £ J dx x dx 2 F^ q , Vi {x x ) F q x ^ q , V2 (x 2 ) a^ v , v ,(x x x 2 s) (67) 

Ai,A2 

This has to be convoluted with the quark structure functions to yield the cross section for the 
pp initial state. 

Eq. (1671) contains integrations over x 12 , the energy fractions of the vector bosons that are 
radiated from the initial quarks, and a sum over vector-boson helicities. In contrast to the 
analogous Weizsacker- Williams approximation for photons, there is a longitudinal polarization 
direction in addition to the two transversal polarization directions. Explicitly, the structure 
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functions are 
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with x = 1 — x. The vector and axial couplings for a fermion branching into a W are 
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The kinematical 



where £3 = ±| is the fermion isospin, and q its charge. 

These structure functions depend on a transverse-momentum cutoff pj_,i 
limit for the cutoff is 

< xv^/2. (71) 

In the derivation of fl68aH68c"I) . one integrates over p± under the assumption that it is small 
compared to the subprocess energy, so the subprocess cross section does not depend on it. For 
the longitudinal structure function that we are most interested in, this can be justified because 
the limit p±, m ax —> 00 is finite. This structure function is concentrated near p± = xmy. The 
transverse structure functions have a logarithmic divergence in p_i_ imax , so the cutoff is needed 
there. This already suggests that the EWA is more reliable for longitudinal than for transversal 
vector bosons. 

In Fig. EH we display the structure functions of W and Z bosons, separately for positive, 
longitudinal, and negative helicity. The emitting quark has been chosen to be an up-type 
quark; for down-type quarks or electrons the Z curves have to be renormalized according to 
the respective charges. For antiquarks or positrons, the transverse polarizations have to be 
interchanged. The plots illustrate the fact that emission of a If or Z, in particular at high 
energies, is more likely for a transversally polarized vector boson. In effect, the production of 
longitudinally polarized VV pairs which couple to the symmetry-breaking sector is suppressed 
compared to this irreducible background. 

Fig. [10] exemplifies the differences between the exact result for qq —>■ qq+VV processes which 
contain resonant weak-boson scattering. To make a meaningful comparison, we first recall that 
in the EWA the initial vector bosons are on-shell, while in the exact process they are off shell. 
The on-shell amplitudes have a Coulomb singularity due to photon and Z, W exchange. In 
particular, an on-shell cross section with photon exchange is infinite, while Z/W exchange 
yields a Coulomb peak proportional to s 2 /M v . Here, s is the cm. energy of the vector-boson 
subsystem, equal to the invariant mass squared M vv of the outgoing vector bosons. To reduce 
this effect which in the exact result is regulated by the vector-boson virtuality, we cut the pt 
of the outgoing vector bosons at 30 GeV. 
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Figure 9: EWA structure functions for W (left) and Z (right) emission from an up-type quark, 
for yfs = 2 TeV and p±, ma x given by ( [71] ). 

A particular choice of this cut allows us to approximate the high-energy end of the Myy 
distribution for the SM with a heavy Higgs (Fig. [!(]] top) quite well [25] . This is misleading, 
however: with the same cut, the prediction of the tensor resonance case (Fig. [TO], lower left) with 
its unitarity saturation beyond the peak is considerably worse. If we are looking at ZZ — > WW 
instead of WW — > ZZ, the EWA background undershoots the exact value (Fig. [TO], lower 
right). More importantly, while the peak can be approximated up to better than a factor 2, the 
background is predicted with less accuracy. Since Myy cannot be reconstructed experimentally 
(apart from ZZ final states), so sideband subtraction is not possible, this significantly affects 
the analysis. 

Part of the deviation is due to the kinematical simplifications inherent in the derivation, 
which can be improved in principle [l9lfT8] . Unfortunately, this only marginally improves the 
EWA, since the main error comes from the existence of irreducible background diagrams for 
on-shell vector boson pair + jets production, and additional irreducible background for the 
complete six-fermion process, cf. Fig. [11] which cannot be accounted for in this way. Off- 
shell, those background diagrams are connected to the signal diagrams by gauge invariance and 
cannot be neglected: simply omitting them would disrupt detailed cancellations, similar to the 
familiar s/t-channel cancellation in W pair production [oTj] . 

5.3 Complete Simulation 

The implementation of the off-shell continued amplitudes in the Monte-Carlo generator WHIZARD 
allows us to get rid of the EWA and to simulate event samples for the complete process pp — > 
qq'+4:f, where the four additional fermions are the decay products of the vector bosons, or come 
from the irreducible background. Using, e.g., PYTHIA for parton showering and fragmentation, 
this results in physical LHC events that can be analyzed by detector simulation and eventually 
compared to real data. 

For illustration, in Fig. [T2] we present the result of a parton-level simulation of WW/ZZ 
scattering, using complete six-fermion matrix elements. In these plots, we compare the effect 
of a 850 GeV vector resonance with the nonresonant (unitarized) LET model, which serves 
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Figure 10: Comparison of the exact (red) and EWA (green, dashed) results for weak-boson 
scattering for processes of the type q\q 2 — > q'iq' 2 VV for \fs qiq2 = 2 TeV. Upper line: scalar 
isosinglet resonance, lower line: tensor isosinglet resonance. The resonance masses and cou- 
plings are Mr = 1 TeV and gn = 1, respectively, the amplitudes are unitarized by the K-matrix 
scheme of Sec. \4-5\ and a pr cut of 30 GeV has been applied to the vector bosons. 



q q q f q q 




q f q q q q 



Figure 11: Feynman graphs that contribute an irreducible background to weak-boson scattering 
in 2 — > 6 fermion processes. E.g. double final state and double initial state radiation, as well 
as t-channel like diagrams. 



29 



as a reference model for the higgsless case. In the four-lepton invariant mass, the resonance 
is clearly visible. However, this quantity is not an observable. The azimuthal distance of the 
two decay leptons is observable; there, vector-resonance exchange in s- and t-channel leads to 
a significant excess. 

A realistic study would be based on a sum over all possible final states with parton shower 
and hadronization, using cuts and distributions in observable quantities. Furthermore, it would 
include a complete account of background and detector effects. A cut-based analysis strategy 
was proposed in Refs. [8|9] . An ATLAS study that makes use of the parameterization of the 
present paper is currently under way (51]. 

6 Summary and Conclusions 

We have described a generic approach to extrapolating vector-boson scattering into the en- 
ergy range where no perturbative predictions exist. Nontrivial features of the amplitudes are 
possible, which will likely appear as resonances. In addition to the classical alternative of a 
heavy scalar-isoscalar (Higgs) or a vector-isovector (technirho or W') resonance, we account for 
scalar-isotensor resonances which are present in extended models, and for tensor resonances that 
could, for instance, be associated with gravity in extra dimensions. Furthermore, we connect 
the model-dependent part to the model-independent low-energy effective theory and keep this 
relation transparent in the implementation. Unitarization of the on-shell amplitudes avoids the 
problem of unphysical behavior at the highest energies that plagues a naive tree-level approach. 

Our approach is economical in the number of free parameters, but intended as a sufficiently 
general description of those energy regions where the LHC will have sensitivity. If necessary, 
refinements of the models, such as recurring resonances or more exotic behavior of the am- 
plitudes, are straightforward to add. The resulting amplitudes are translated into effective 
form-factors for vector-boson vertices in unitarity gauge. This allows for an implementation 
in universal Monte-Carlo event generators, which we have realized for the case of the WHIZARD 
event generator. 

While the leading electroweak loop corrections for vector-boson scattering are included, 
QCD corrections are not yet implemented. These have been considered in Ref. |52j and should 
be combined with the effects modeled by our approach. 

With the event generator at hand, model-independent studies and analyses of vector-boson 
scattering, both in SM extensions and in Higgsless models, become feasible. No approximations 
beyond those inherent in the modeling are involved, as it is essential for unbiased data analysis. 
A particular feature of our implementation is the smooth transition to the SM case (with a 
Higgs boson) or, alternatively, to a featureless LET model of strong WW scattering without 
resonances, respecting unitarity. In data analysis, the signal can be defined as the deviation 
with respect to either one of those reference models. 
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Figure 12: Unweighted event samples for pp — *> ud + e + v e e~u e . Upper: Invariant mass of 
the e + v e e~v e system. Lower: Azimuthal distance of the charged leptons. The red histogram 
(hatched) corresponds to a model with a vector resonance (M p = 850 GeV and g p = 1). The 
green histogram (filled) is the LET result. Both models have been unitarized by the K-matrix 
scheme. Cuts: p\_{tv) > 30 GeV, |?7(^)| < 1.5, 6{u/d) > 0.5°. The integrated luminosity is 
4 x 225 fbT 1 (the factor 4 accounts for the sum over e, fx). 
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A Conventions and algebra 

A.l 577(2) algebra 

Throughout this paper, we use boldface notation for objects that are defined in the adjoint of 
SU(2), e.g., 

W M = Wf-^ (72) 

with the Pauli matrices r a (a = 1,2,3), and summation over a understood. 

For describing isospin quintet resonances, we introduce tensor products of Pauli matrices: 

r ++ = r + ® r + (73) 

r + = l( r + ® T z + T z ® r +) (74) 

r° = ^g(r 3 <g>T 3 -T + <g>T~ -t~ <g>r+) (75) 

t~ = \{t~ ®t 3 + t 3 ®t~) (76) 

r~~ = r~ ® r~ (77) 



These are normalized: 
Isospin singlet: 



tr [r ++ r— ] = tr [r+r"] = tr [r°r°] = 1 (78) 



r aa = r a ®r a = r 3 ® r 3 + 2r+ <g> r~ + 2t _ <g> r+ (79) 
Tracing this with something else gives 

tr [(A <g> 5)r aa ] = 2 tr [AS] , (80) 

in particular 

tr [r ++ r aa ] = tr [r + T aa ] = tr [rV a ] = 0, tr [ r 33 r aa ] = 4, tr [r aa r bb ] = 12. (81) 
Furthermore we need: 

tr [r°(r 3 ® r 3 )] = -i= tr [r°(r+ ® r" + r" ® r+)] = (82) 

A. 2 Goldstone bosons and Gauge fields 

We define the Goldstone scalar triplet w 1 ' 2 ' 3 or, alternatively, w + ,w~,z such that 

w 1 = -^(yW + + w~) w + = -^(w 1 — i^ 2 ) 

w 2 = -L=( w + - w ~) w - = + i w 2 ) (83) 

and w 3 = z. Contractions: 

w = w a r a = V2(w + t + + W'T') + zt 3 (84) 
(w) 2 = w + w~ + w~w + + zz (85) 



32 



The Higgs-field matrix is given by 

£ = exp 

The covariant derivative of the Higgs field is 



-w 



(86) 



DT, = <9£ + i 3 WS - u/E B 



(87) 

Unitary gauge would mean w = 0, i.e., E = 1. Herewith, we define the vector field 

V = £(£>£)t = -(DS)S t , (88) 

which is in the adjoint representation of SU(2) L , and is a linear combination of Pauli matrices. 
Hence, tr [V] = 0. Note that V is antihermitian, V" 1 " = —V. 

Gauge fields for the electroweak and strong interactions are defined such that they transform 
under SU(2) L x U(1) Y as W -> U L WU[, 

= d,W„ - d„W„ + ig[W„, W„], (89) 

= E (d^B u - d u B^) ySt; (90) 

furthermore there is the QCD gauge field = G^J^. 

In the gaugeless limit, the expansion in terms of Goldstone fields is V =>- ^dw k r k + 0(v ~ 2 ). 
Expressing this in terms of charge eigenstates, we derive 



V 



1 



V2dw + r + + V2dw- T - + dzr 3 
V2w + dzr + — \[2w~ dzr~ — w + dw~~T 3 



+ 0{v~ 3 ) 



and thus 



Hence, 



tr [VpV,,] = -— (d^w + d u w + d^w d v w + + d^zd v z) + 0(v" 



1 

tr [V • VI = dw + dw + -dzdz 

4 2 

In the notation used for couplings to isospin quintets, we have 
^V {M <g> V,} = - ^2d f ,w + d u w + r ++ + 2d lt w-d v w-T— 

+ \/2 (d fl w + d u z + d u w + d fl z) r + + a/2 [d jJ w~d v z + d u w~d^z) t~ 



(91) 



(92) 



(93) 



+ d^zd v zr 3 <g> r 3 + (d^w+dvW + <9^u> (9^w + ) (r + ® r + r <g> r + ) 



(94) 



And, 



41 2<9w + • <9w + r ++ + 2duT • dw~T~ + 2v / 2<9w + • dzr + + 2 v / 29w" • 9zr" 
v 2 { 

+ dz ■ dzr 33 + 2dw + ■ dw~ (t + ~ + r~ + ) j (95) 
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A. 3 Tensor Fields 

A massive tensor field f^ v is subject to the conditions 

r = r, /% = o o,r = d„r = o. (96) 

Its spin sum is given by 

t *^ v t P° — 1 ( pw p vu _|_ pP a p v P^ — - ( pw pP a } (97) 
a A A 2 3 

where 

hPk u 

pi»{k) = gr-- W . (98) 

The free Lagrangian is 

M 2 

C, = C kin - —f^r (99) 
where the kinetic part corresponds to the spin sum (f97|) . 

A. 4 Integrals in spin-isospin eigenamplitudes 

To get compact expressions for the spin-isospin eigenamplitudes, we define the following inte- 
grals: 



Sj{s) = j jj-^P (t,s,u)Pj(s,t,u) (100a) 
f° dt t 

Vj{s) = J jj-^P 1 {t,s,u)Pj{s,t,u) (100b) 

r° dt t 2 

Vj{s) = J -—^P 2 (t,s,u)Pj(s,t,u) (100c) 
The integrals over u 2 /(u — M 2 ) are (— 1) J times those over t 2 /(t — M 2 ). 
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Explicit expressions for these integrals are: 
<s A/f 4 M 2 

*= M S t 7 bt rf (101a) 

M 4 <s M 4 M 2 
Si(s)= 2— + U— (2M 2 + S )log— — (101b) 

s 6 s 2 s + M 2 

M 4 /W 4 M 2 

S 2 (s) = -y (6M 2 + 3s)+ — (6M 4 + 6M 2 s + s 2 ) log -— — (101c) 

s + M 2 

M 4 M 4 M 2 

S 3 (s) = — (60M 4 + 60M 2 s + lis 2 ) + -^(2M 2 + s) (lOM 4 + 10M 2 s + s 2 ) log — 

3s d v s 4 s + M 2 

(lOld) 

9 S + M 2 M 2 

s s + 

M 2 -4- 2 <? / M 2 \ 

^(5) = — 1— (2s + (2M 2 + s) log -^j (lOlf ) 

1 1 M 2 

V (s)= -(2M 2 + lls) + -(M* + QM 2 s + 6s 2 )\og—— (lOlg) 

2 s s + M 2 

i r m 2 

^(s) = _ <^ S (12M 4 + 72M 2 s + 73s 2 ) + 6(2M 2 + s)(M 4 + 6M 2 s + 6s 2 ) log 

6s 2 [ s + M 2 



(lOlh) 



B Feynman rules for scalar and tensor resonances 



We briefly summarize the Feynman rules for scalar and tensor resonances that derive from the 
interaction Lagrangians (llOal) - (llOep . The ks in this section are the momenta of the Goldstone 
bosons. 

Scalar isoscalar: 

, _ 2ig 2z<7, 
aw w : (k + ■ k-) azz : ik\ ■ k 2 ) (102) 

v v 

Scalar isotensor: 

: - ^-(h ■ k 2 ) (jPvFz : - -(fc, ■ k z ) (103) 

v v 

<fP zz : - ^L(k! ■ k 2 ) <f) w + w- : + • fc_) (104) 

V 3w V 3f 

For the Feynman rules of the tensor resonances we use the symbol C^ pa := g^pgua+g^gup — 
\gp V g P a to get (momenta incoming). 



Tensor isoscalar: 



f^ u w + w : —C 'au Mok+k ^ f^zz : -C 'n U po -^i^2 (105) 

v v 

Tensor isotensor 



a^wV : - ^C^K a%w^z : - % ^C^k%K (106) 

a pv zz '■ ~7F~Cp,v,pok' f \ k 2 a^w^w : + - y= C jlvpa k p A k_ (107) 
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Note that taking the conditions on the tracelessness as well as the transversality not nec- 
essarily demands the coupling of the tensor resonance to a conserved current (like the energy- 
momentum tensor) which leads to the same Feynman rules as in [53]. The constraint of the 
LET on the other hand results in an (off-shell continued) amplitude that is identical to the one 
of a massive graviton resonance. 

C Vector Resonance Exchange 

Heavy vector resonances have been studied many times in the literature, and various different 
formalisms describe their interactions with the SM particles. In this section, we demonstrate 
the equivalence of some popular approaches. In particular, we look at the correction to the 
amplitude A(s,t,u) for Goldstone-Goldstone scattering which via the GBET and spin/isospin 
symmetry yields the leading term for all channels of quasi-elastic WW scattering, w + w~ — ► zz. 
Since we maintain manifest SU(2)l x U(1)y gauge invariance by using covariant derivatives, 
the GBET holds in any formalism that we describe. If desired, this can be verified by switching 
to unitarity gauge and computing the W/Z scattering amplitudes directly. 

We ignore couplings to fermion currents, which generically will be present and sizable. 
While such couplings get shifted by some of the transformations described below, in our model- 
independent approach they are considered as independent parameters which are determined 
by independent measurements. The shifts of fermion couplings induced by reparameterizations 
merely change a set of undetermined parameters into another set of undetermined parameters, 
so for our purposes there is no need to calculate them. However, in the context of a specific 
model, one should always treat fermionic and bosonic sectors together when applying reparam- 
eterizations [54] . A specific example for this in the context of Little Higgs models can be found 
in [55] . 

1. We use the representation with the Goldstone field S = exp(— iw/v) and introduce the p 
resonance as a vector field in the iso-triplet representation. The Goldstone kinetic term 



is 





and can be expanded up to second order as 



(0) 



dw + dw H — dz dz, 
2 

(w + dw~ + w~dw + + z 

2v A 



(109) 



kin 



kin 



dz) 2 . 




The interaction Lagrangian is 



2 



,2 



int 




(111) 
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Figure 13: Feynman graphs for w + w~ — > 2:2 : contact term from the non-linear Lagrangian 
leading to the LET on the left, resonance t- and u-channel exchange in the middle, and contact 
terms from p/w/z mixing on the right. 

Ignoring gauge fields, we carry out the expansion up to third order: 
CQ = -g p v (p + dw- + p-dw + + p°dz) , (112) 
Ant = ^9p [p + ( w ~9z — zdw~) — p~ (w + dz — z dw + ) + p° (w + dw~ — w~dw + )~\ , (113) 
= [p + ((w + w~ + z 2 ) dw~ — w~ 2 dw + — w~ z dz) 

+ p~ ((u> + u>~ + z 2 ) dw + — w +2 dw~ — w + zdz) 

+ p° (—w~ z dw + — w + zdw~ + 2w + w~dz)~\ , (114) 

leading to the corresponding Feynman rules 

p^w^ : igvk T4l p°w + i/T : igvk z ^ (115) 

p p w T z : ± ig(k z - k T )^ pJtu+uT : ig(k_ - k + ) p (116) 

p^w T zz: - i-^-{2k T - k x - k 2 ) tl p°w + w~z: - i^-(2k z - k + - k^) p (117) 

The resulting Feynman graphs for w + w~ — > zz are shown in Fig. [131 There is a Goldstone 
contact interaction from the Goldstone kinetic term; this yields the low-energy theorem 
(LET) A(s,t, u) = s/v 2 . Resonance exchange adds t- and it-channel contributions, 



u s — t 
t-M 2 ~ + u-M 2 



A' p (s, t, u) = -g 2 p ( + —1, ) , (Hi 



which in the limit s — > become A' p — > —3g 2 s/M 2 . Furthermore, there are contributions 
where the resonance mixes with the external Goldstone scalar, either as an external 
wave-function correction, or with the p coupling to three scalars. In both graph types 
the resonance pole cancels, so they are additional contact terms proportional to s/M 2 . 
In unitarity gauge, this translates into a correction to W/Z masses and interactions. 

The appearance of a contact term looks like a violation of the LET. However, the measured 
value of v (e.g., as extracted from W/Z pole data) includes those additional contributions, 
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so they merely renormalize a fictitious bare v value. This renormalization can be made 
explicit by adding a counterterm to the p interaction Lagrangian, which by power counting 

2 

and symmetries must be of the form af^C^n with an appropriate prefactor a. In effect, 
expressed in terms of the observed scale v, the LET holds, and the vector-exchange 
amplitude is given by 

(s — u s — t s \ 

T^W + ^W + 3 Tp)> (119) 

which vanishes as s 2 as s — > 0. 

2. In the previous paragraph, the vector resonance was coupled to W/Z bosons by a mass 
mixing term, tr [Vp] . Alternatively, we could couple it by a kinetic mixing term, 

£int = -2ptr[W^ P n (120) 

where the resonance "field strength" is p^ v = D^p v — D u p fl with the covariant derivative 
in the adjoint representation. Partial integration gives 

£ int = 2^tr[(^W^)pH- (121) 
Here, we can apply the W field equation 

D»W^ = ijV v (122) 

to obtain 

Ant = ^tr[V^] (123) 

as before, so we get the same scattering amplitude. Using the equations of motion is 
precisely an application of the UET. 

3. In the CCWZ formalism [31], the elementary building block is £ with = S. From £, 
we can construct a vector and an axial vector field, 

v M = \ {?d^ + e^e f ) and a> = ~ (e f £u - e^e f ) • (124) 

Under SU(2) C , these transform like a gauge field and a matter field, respectively, 

V ^ U c VU j c - (D^UcWl and A^U C AU ] C . (125) 

A is related to the vector current that we have used in our previous formulation: £'V^£ = 
2iA^. We just have to redefine p M — > to obtain 

£ kin = -2t; 2 tr[^^] and £ int = tr [p^] , (126) 

so a matter field p coupled to the axial vector A yields the same scattering amplitude 
again. 
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4. Alternatively, we can couple p to the vector field V by assigning to it a gauge- field 
transformation law under SU(2)c, 

p -> UcpUl - i 2 -ML{D,U c )Ul (127) 
so the leading invariant term containing p is 

M 2 

= -g 2 p v 2 tr [VV] - ig p vM tr [Vp] + — tr [pp] (128) 

The expansion of V is even in the number of Goldstone fields. Therefore, in this expression, 
the last term is the p mass, the second term yields the p°w + w~ and p ±r w T z couplings, 
and the first term is a contact term. Note that the gauge coupling is proportional to 
l/g p - The resulting w + w~ — > zz amplitude is again (11191) . this time without the need for 
renormalizing v. 

5. The BESS model [6] has a similar setup. Instead of gauging just SU(2)c, we can extend 
the local symmetry by an extra local, nonlinearly realized SU(2)l x SU(2)r. This results 
in two vector isotriplets, which can be combined to a vector and an axial vector isotriplet, 
respectively. Only the vector couples to longitudinal W/Z pairs, and the amplitude (11191) 
can be derived as before. 

The different formalisms for coupling vector resonances all result in the same scattering ampli- 
tude. This is not surprising since this amplitude is completely determined by spin and isospin 
conservation together with the LET. In order to give the CCWZ interpretation of the vec- 
tor resonance as a gauge field (in contrast to a generic matter field) a physical meaning, we 
would have to measure triple p couplings, analogous to the LEP2 measurements of triple gauge 
couplings. Unfortunately, such measurements are beyond the reach of LHC. 

D Specific Models 

In the literature, a variety of "benchmark" models has been formulated that test weak-boson 
scattering. In this section, we relate some of them to our parameterization: 

1. The SM. As discussed in the main text, for g a = 1 the scalar resonance model pre- 
cisely reproduces SM Higgs exchange. Alternatively, one can switch to the default SM 
implementation (in WHIZARD) without extra resonances. 

2. Scalar resonances. Refs. [8.9J introduce a collection of models, among them two with a 
scalar resonance ("0(2 AT)" and "chirally coupled scalar"). The latter model is identical 
to our scalar resonance parameterization. The 0(2N) model is essentially a special case 
of this with fixed mass and width; the only distinction is a logarithmic cutoff-dependent 
modification, which manifests itself beyond the resonance. This detail is unlikely to be 
detectable at the LHC. 

3. Vector resonances. The chirally-coupled vector resonance model of [8f9] is identi- 
cal to ours (see the discussion of the CCWZ formalism in App. where we identify 
a = {2g p v / M p ) 2 and g = M 2 /{2v 2 g p ). An analogous identification holds for the BESS 
model [511, with a replaced by f3 in their notation. 



C 



hit 



2 2. 

-g p v tr 



V 



M 
i p 

2g P v 
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Pade/IAM unitarization model. 

As discussed in Sec. I4.5[ this scheme is a special case of the K-matrix scheme as defined 
in the present paper. For a given combination (a 4 , a 5 ) we use Eqs. (I38al I38c[ I38el) to 
determine the NLO correction Ayj to the three amplitudes Aqq, An, and A20 which 
without correction would violate unitarity. Then, we can use ( )56|) to identify scalar, 
vector, and tensor resonance masses and widths. If we neglect the loop corrections in 
f l38aH38e j ) . we obtain 

Ml = r CT = ^ (129a) 

4(7a 4 (/x) + Ha 5 (/i) 16tt y j 

M " = 4(a 4 (/x) - 2aM) Tp= 9^ (129b) 



Ml = , ^ „ T f = (129c) 

J 16(2a 4 (» + a 5 (ji)) f 32n 



where we have to define a renormalization scheme and fix the scale /1. Note that the 
tensor-resonance parameters are unphysical. This is due to the negative sign of A$ in 
Eq. (1371) . This model ignores the possibility of isotensor resonances <p or a. 



E On-shell vector boson scattering 

In this section we summarize the plots for "partonic" scattering of spin-averaged and summed 
vector bosons. In all these pictures, the EW gauge bosons are treated on-shell, hence the cross 
sections start when the physical WW or ZZ threshold is reached. Since we did not switch 
off the electromagnetic coupling in those plots, we applied a cut of 15 degrees around the 
beam axis to cut out the Coulomb scattering part. Fig. [TH shows in the upper line the SM 
with a 120 GeV Higgs on the left and a heavy 1 TeV on the right. Unitarity is preserved in 
those cases because of the (s-channel) Higgs exchange. Besides the dominant resonance for 
a heavy Higgs, the amplitudes show a saturation for the high-energy tails which starts again 
violating partial-wave unitarity for 1.2, 3.5, and 1.7 TeV for the I = 0,1,2 isospin channels, 
respectively. Completely removing the Higgs as in the middle line of [TH leads to a rise of the 
amplitudes (and hence the cross sections) with s (the z z — > zz process is absent in that case). 
Switching on the A-matrix unitarization damps the amplitudes back to the Argand circle, 
thereby restoring unitarity. This happens for the above mentioned values for the corresponding 
isospin eigenamplitudes. In the lower line of Fig. [TH the case of the LET extended by nonzero 
values for the parameters a 4 and a*, are shown, on the left the badly diverging case without 
unitarization, and the K-matrix unitarized case on the right. 

In Fig. [Tj)]we show the cross sections for the five different vector boson scattering processes 
with the presence of the five isospin-allowed resonances mentioned in the text. ZZ — > ZZ and 
WW — > ZZ contain all three isospin eigenamplitudes, hence show a resonance in all channels 
except for the vector resonance case where it is forbidden by the Yang-Landau theorem. The 
amplitude WZ — > WZ does not have isoscalar resonances, while the amplitude W + W + — > 
W + W + allows only isotensor resonances. Finally, W + W~ — > W + W~ contains all resonances. 
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Figure 14: Cross sections (in nanobarns) for the five different scattering processes of longi- 
tudinal weak gauge bosons: SM with a 120 GeV and a 1 TeV Higgs in the upper line, in the 
middle: SM without a Higgs without and with K-matrix unitarization, respectively. In the lower 
line, the case of switched on are shown, on the left without, on the right with K matrix 
unitarization. The contribution from the forward region is cut out by a 15 degree cut around 
the beam axis. 
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Figure 15: Cross sections for VV — > VV scattering in nanobarns, with the presence of reso- 
nances: scalars (is o scalar a and isotensor <f>) in the upper line, vector isovector p in the middle, 
and tensors (isoscalar f and isotensor a) in the lower line, respectively. All amplitudes have 
been unitarized according to the K-matrix method. The resonance mass is set to 500 GeV in 
each case. Again, the contribution from the forward region is cut out by a 15 degree cut around 
the beam axis. 
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